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On the Shapes of the Isogeotherms under Mountain Ranges in 

Radio-active Districts, 
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(Received January 7, — ^Eead February 10, 1910.) 

The importance attached by geologists to the distribution of temperature 
within the earth's crust as a factor in the production of movements of the 
crust, and in particular in the formation of mountain ranges,* has made it 
necessary to consider if it is possible to determine the distribution of 
temperature under and in the neighbourhood of a mountain range, by a 
method more rigid and accurate than that used by Fisher,*]- and more closely 
following the physical condition of the problem than that used by Thoma.J 

In what follows it will be shown that an accurate solution can be obtained 
in certain simple cases, even when the soil is radio-active. 

The problem as presented in nature is too complicated to admit of mathe- 
matical treatment, and to render it more manageable the following 
^simplifications are introduced : — 1. The curvature of the surface of the earth 
is neglected. 2. The mountain range is assumed to be of uniform height 
and shape, and to have a length great in comparison with its height and 
width. 3. The section of the range is taken to be one of a few simple and 
typical forms. 4. The material of the range and that of the earth's crust 
in the neighbourhood are taken to have the same radio-activity and the 
same conductivity for heat.§ 5. The isothermal surfaces (isogeotherms) 
where there* are no mountains are assumed to be planes parallel to the 
surface of the earth. 

Taking the axis of x horizontal and at right angles to the axis of the 
mountain range, and the axis of z vertically downwards from the surface of 
the plain at the foot of the range, we have the following equations which 
must be satisfied by the temperature v at the point ^ ^ of the crust, 

S'4? + 9lVH=:0 (1) 

and 

/.g = M.-Y) (2) 

^ See, e.g., Prof. Joly's ^ Eadio -activity and Geology,' chap. 5, p. 101. 
t Fisher, ' Physics of the Earth's Crust,' p. 222, 2nd edit. 
J Thoma, ' Dissertation,' Freiburg, 1906 ; ' Beiblatter,' vol. 31, p. 850, 1907. 
§ Sandstone, slate, limestone, marble, dolomite and granite all have conductivities 
•about 0*005. Basalt has a higher, syenite and chalk have lower conductivities. 
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at the surface, where h is the heat conductivity of the material of the 
crust; ID is the heat generated per cubic centimetre per second by radio- 
active processes ; n is measured inward along the normal to the surface ; 
V is the temperature of the air in contact with the surface at the point ^ 
and h is the external conductivity or emissivity of the surface. 

The complete integral of equation (1) consists of the complementary 
function, i.e. the solution of 

|| + |l = (l'> 

and a particular integral of the original equation, which in this case ia 
obviously —w/k.^z^. 

The variation of the temperature of the free air with altitude may be taken 
as uniform throughout the lower portion of the atmosphere, and of the order 
0°*5 C. per 100 metres of ascent, and it may be shown that for surfaces of 
which the angle of slope is not great the temperature of the surface will vary 
in the same simple way. Actual observations of temperature on mountain 
slopes confirm this deduction, and show that there is a decrease of tempera- 
ture of between 0^*54 and 0^*66 0. per 100 metres of ascent.* 

For the present purpose it will therefore be sufficient to substitute for 
equation (2) the following, 

dv/dz = a' (2'> 

at the surface of the mountain, where a is a constant. 

For a given contour of section of the mountain range, the above equations 
provide sufficient data for the solution of the problem, but there are still 
difficulties of a mathematical nature in obtaining the solution. These diffi- 
culties are obviated if, instead of determining the temperature for a given 
shape of section, we write down the expressions for the temperature which : 
(a) satisfy equation (1), (h) are symmetrical about the axis of z, (c) have 
maximum values at a? = 0, (6^) become zero at x or z = infinity, (e) have no- 
singularities within the region considered, and determine the shape of section 
for which these expressions or any combination of them form the accurate- 
solution. 

The expressions which satisfy these conditions are 

A^ £- log {x^ + (^ + a^y} , 

where n is an integer. An and a^ constants. 

To these must be added the further solution v = az, where cc is a constant,, 
and the particular integral. 

■^ Hann, ' Handbucli der Meteorologie,' pp. 101 — 104. 
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Hence the temperature throughout the section is given by 

^-Vo = az-f . iV" 2°°A„^log {x^ + {z + a„f}, (3> 



A/ ju n — n C(y^ 



n = 



where Vq is the temperature of the surface of the plain at the foot of the 
range. 

The temperature of the surface of the mountain at different levels is found 
by integrating equation (2'), and is given by the equation 

t;— Vo = oc'z, (4)' 

Since at the surface of the mountain both equations (3) and (4) hold, the 
equation to the contour of the range is obtained by eliminating v from the 
two equations. We thus obtain 



= (^_^')^_| . 1% I A,|llog {x'-]-(z + anyy^ (5> 

/^ ^ •j^ = ^^ 



as the equation to the contour of the section of the range. 

By a proper choice of the constants of the types An and an, any contour 
symmetrical about the axis of z, and not having singular points, may be 
represented by equation (5), then equation (3) gives the distribution of 
temperature throughout the section. 

The determination of these constants for a particular form of contour is, 
however, not a very easy problem, and it is doubtful whether it is worth the 
trouble of representing the contour very accurately when the material of the 
section deviates so much from the homogeneous medium assumed in 4, p. 339.. 

A considerable amount of useful information may, however, be obtained 
from the theory given, by taking so simple a form of contour that the series 
of differentiated terms reduces to a single term. The term which gives a 
contour most closely resembling that of a mountain range is the one for 
which 7^ = 1, which gives as the equation of the contour — 

= (cc-cc')z-f . $ + A~ ^ + "^ ., (5'> 

^ "^ k 2 x^ + (z + af ^ ^ 

and for the isooeotherms — 



2 



■XT w z'^ ^ . z-\-a ,o/\ 

k 2 x^-^(z-\-af ^ ^ 

^ For a series of parallel ranges with peaks at ^j, Xo,^ etc., the equations become 

W ^ CO 00 ^n 

'^ ^ m = lw:;=0 ^^ 

?/' '~^ oo 00 rj"^ 

and = (a-a')^--7.^+ 2 2 A,^«;,-nlog {(^-^^0' + (^ + a„,>}.. 
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In the equation to the contour write -^/(oc — u) = /3 and A/(oi--a) = B, it 

then becomes 

z (l-/3z) + B -^P^:^ = 0. (6) 

At ^ = — H let ^ = 0, and let H be numerically the smallest value of z for 
which this is possible ; that is, let H be the height of the range. Then 

-H(l+/3H) + B/(a-H) = 0. (7) 

Further, let ^ = & when z = —h, wiiere h lies between and H. Then 

TT /t (1 + ^h) ^ (c(,-B.)(a-h) _ w'-('K + h)a + B.h 

"^^ H(1+;SH) b^ + {a-hy d'-2ha + ¥ + ¥ ' 

^^^ H(1+/SH)-^- ^^^ 

The equation becomes 

a2_2Ha + 2H/i-(/^H&')== 0, 

which gives a == 'R-{-^{{ll-^hy-\-'b^), (9) 

the negative value of the radical being discarded, as it would give a singular 
point within the contour. 

Finally, we have B - H(1 + /3H) v'{(H-/0' + &2}„ (10) 

Hence the equation to the contour may be written down if we know the 
height H of the range and the half breadth 5 at a height h, nearly half that 
of the range. 

If we write d for the distance from the summit of the section of the range 
to the point on the surface whose height li above the plain is determined by 
equation (8), then 

and we get 

^=: H + 6Z, (90; B = H(l + /3H)f^, (10'); A = (a-»^0 H^(1 + /3H). (11) 

Hence, if a long mountain range of uniform height H and the plains 
at its base be composed of material of heat conductivity h and of radio- 
activity such that 10 gramme-degrees of heat are generated in each cubic 
centimetre per second, and if ol be the increase of temperature per 
centimetre of descent into the material of the plain, and a' the decrease 
of temperature per centimetre of vertical ascent up the sides of the range, 
and if the height h of a point on the surface of the range be determined so 
that if /3 be written for iol2h {oc-a), then h (l + ^h)/!! (1 + ^H) = |, and the 
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length of tlie line joining this point to the summit of the section be d, the 
equation to the contour of the section may be written 

z(l^l3z) + (a--a')lid(l + 0m- ^^ ^ + ^ + ^ ... = 0, (5") 

and the equation to the isogeotherms 

.-Yo = ..-I . |V(.-.0H.(l + /3H)^,;ff+^^^, (3") 

Subject to the conditions that ^ be not less than H (l + ^H)/(l + 2/3H),* 
-and that z be not taken greater than 1/20, z being measured downward from 
the surface of the plain, and x measured from the centre of the section of 
the range. 

The temperature within a horizontal tunnel through the range at a height 
Z above the plain is given by the equation 

^,_To = ---aZ-'^|%(a~aOH^(l + ^H)-— -5±^^ . (12) 

k 2 ^ / V -TA^ '^^2 + (H + r^-Z)^ ^ ^ 

The excess of temperature at a point at the centre of the section over that 
at a point on the surface at the same height is 

(^^^')^_^ . |' + (a~-aOHrf(l + ^H)(^+H + ^)"^ (13) 

When the point is level with the plain this excess becomes 

(«-,') (l+;9H)j^^. (13') 

The rate of increase of temperature vertically downwards at the point z 
in a vertical line through the summit is given by the equation 

||=^-f^-(^-^')Hc^(l + /3H)(^ + H + rf)-2, (14) 

which at the summit itself becomes 

|^= a+|H-^(^~-aOH(l+^H)/d (140 

The rock constants which enter into these equations will have different 
values in different districts, but it is worth while working out in detail an 
ideal case in which they have values approximately the mean of the values 
found experimentally, in order to realise the relative magnitudes of the terms 
of the final equations. The mean rate of increase of temperature downwards 
from a level surface may be taken as 1° C. in 33 metres,t which gives 

■^ This is the condition that H is the least root of equation (6). 

t ' Brit. Assoc. Rep.,' 1882, p. 88 ; Prestwich, ' Roy. Soc. Proc.,' vol. 41, p. 1, 1886 ; 
Joly, * Radio-activity, etc.,' p. 71. 
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a = 3*0 X 10"^ degree per centimetre. The rate of decrease of temperature 
up the sides of mountains varies from 0*^*54 to O'^'GG C. per 100 metres of 
ascent;^' hence a may be taken as 6*0 x 10""^ degree per centimetre, and a — of 
as 2'4 X 10~* degree per centimetre. 

The heat conductivity of rock may be taken as on the average 0'004.f 

The estimates of the mean radium content of rocks differ considerably from 
each other. Strutt^J Eve and McIntoshj§ and more recently Farr and 
MoranceJI give figures of the order 1*1 x 10"*^^ gramme of radium per 
gramme of rock, while JolylT gives 5*0 x lO""^^. 

Taking for the present the smaller estimate, and taking the average density 
of the rocks investigated as 2*5, we have 2*7 x 10"^^ gramme of radium 
per cubic centimetre of rock. Taking Eutherford's estimate*'* of 5*6 x 10~^ 
gramme-degree as the heat given out per second by 1 gramme of radium in 
equilibrium with its associated elements, we get 1*5 x 10~^^ gramme-degree 
generated per second in 1 cubic centimetre of rock as the value of w. Hence, 
^ = ta/2J€(a--'a) is equal to 8 x 10""^ and /3H for ranges of 4 kilometres in 
height will have the value 0*03. Thus h may be determined with suificient 
accuracy for ranges not excessively high if equation (8) be written 
h = ^N (1 + J^H), (8'). From this we see that h is 0*4 per cent, greater than 
half the height of the range for each kilometre of height of the range, and 
that the factor (1+ySH) in the equations (3'') and (5") only differs from 1 by 
0*8 per cent, for each kilometre of height of the range. 

Thus a small amount of radio-active matter in the rocks does not seriously 
modify the shapes of the isogeotherms under ranges of moderate height, and 
in the following graphical representations of them in a few typical cases it 
will be neglected. 

In Diagramsft 1, 2, and 3 the values adopted for the constants are : 
^ = 3*0 X 10"^ degree per centimetre, a — 6*0 x 10"^ degree per centimetre, 
temperature of the plain 14° C, H = 2 kilometres, d = H, -Jx/S H, and 
1^/17 H respectively, which correspond to & — ^^/S H, H, and 2 H respectively. 
According to equation (7), H is the least value d can have. For that value 
the slope of the surface of the range is 45*^ at the summit and decreases. 

'^ Harni, ' Handbuch der Meteorology/ pp. 101 — 104. 

t Joljj ' Radio-activity, etc.,' p. 74. 

J Strutt, * Roy. Soc. Proc.,' A, vol. 77, p. 472, and A, vol. 78, p. 150, 1906. 

§ Eve and Mcintosh, ' Phil Mag.,' vol. 14, p. 231, 1907. 

11 Parr and Florance, ' Phil. Mag.,' vol. 18, p. 812, 1909. 

IT Joly, ' Biit. Assoc. Rep.,' 1908, p. 683. 
^^' Quoted by Joly^ * Radio-activity, etc.,' p. 11. 

ft I have to thank Mr. G. D. West, one of my senior students at the East London 
College, for drawing these curves from the equations. 
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downwards. For greater values of d the slope is less. According to Sir 
Archibald Geikie,* slopes exceeding 40^ seldom exist except in solid rock. 

It will be seen that the isogeotherms under the range are convex upwards^, 
as one would expect from elementary considerations.f The excess of 
temperature at a point under the summit of the range and at the same level 

H 

as the plain over that of the surface of the plain is (a — a) - — -:— ,~ if the 

^ ^ ^ ^ l+H/<i 

rock is not radio-active, and is given by equation (13') if the rock is radio- 
active. The vertical temperature gradient at the summit is 

^(l-H/6^)-f ^'H/^ 

if the rock is not radio-active, and is given by equation (14') if it is radio- 
active. The calculations of Fisher/I which led him to conclude that the 
shape of the surface of the range has little effect on the temperature gradient, 
and that a mountain has solid roots beneath it projecting downwards into 
a substratum of molten rock, are therefore not confirmed by the present 
work. 

^ Geikie, ' Text-book of Geology,' vol. 1, p. 52, 1903. 

t See Everett, ' Brit. Assoc. Rep.,' 1879, p. 345. 

X Fisher, ' Physics of the Earth's Crust,' pp. 222—229, 2nd edit., 1889. 



